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An extension to the Soft-Collinear-Effective Theory (SCET) description of hard jets is motivated 
to include the leading contributions between the propagating partons within the jet with partons 
radiated from a dense extended medium. The resulting effective Lagrangian, containing both a 
leading and a power suppressed (in the hard scale Q 2 ) contribution, arises primarily from interactions 
between the hard collinear modes in the jet with Glauber modes from the medium. In this first 
attempt, the interactions between the hard jet and soft and collinear partonic modes have been 
ignored, in an effort to focus solely on the interactions with the Glauber modes. While the effect of 
such modes on vacuum cross sections are suppressed by powers of the hard scale compared to the 
terms from the SCET Lagrangian, such sub-leading contributions are enhanced by the extent of the 
medium and result in measurable corrections. The veracity of the derived Lagrangian is checked 
by direct comparison with known results from full QCD calculations of two physical observables: 
the transverse momentum broadening of hard jets in dense media and a reanalysis of the transverse 
momentum dependent parton distribution function (TMDPDF). 



I. INTRODUCTION 



The study of hard jets in QCD is now a considerably mature science. Experiments at e + e~ annihilation, Deep- 
Inelastic Scattering (DIS) and p-p machines have yielded a wide array of measurements on a variety of jet observables 
including single particle production, multi-particle correlations as well as event shapes. On the theory side, sophisti- 
cated factorization theorems have been written down which factorize the final state "jet function" from the initial state 
and the hard cross section at leading twist [l[ . In the derivation leading to such factorization theorems 0, H, 0, HI, @| 
the infinite class of Feynman diagrams are subjected to a Landau analysis. Regions of momentum space which yield 
pinch singularities are identified. These represent the leading contributions to such processes and may be decom- 
posed into classes of Feynman diagrams which in turn allow for proofs of factorization. An alternate and equivalent 
approach has recently been afforded by the methods of effective field theories such as the Soft-Collincar-Effective- 
Theory (SCET) 0,1, & dUUH. While not specifically devised to re-derive factorization, SCET presents a formalism 
where the analysis resulting in the identification of the leading contributions may be carried out within the QCD 
Lagrangian resulting in the derivation of an effective Lagrangian which is only applicable to processes within the 
prescribed kinematic regime. In such a formalism, factorization occurs at the level of the Lagrangian and at the level 
of operators fl2l . [l3j . The Feynman rules which arise from an expansion in a small parameter A may then be used to 
systematically study hard processes. 

Power corrections to hard processes in vacuum are suppressed in the presence of a hard scale Q 2 3> ^qcd- However, 
there exist scenarios where a specific set of power corrections (often arising from operators with higher twist) may be 
enhanced and become non-negligible compared to the leading process. One example of this is the case of single spin 
asymmetry in DIS, where leading twist processes yield vanishing results [lil ]. Another example, where the leading 
twist term does not vanish but power corrections may become enhanced is in the presence of a medium is in the case 
of DIS on a large nucleus [l5|, [l6| . The inclusive cross section receives contributions from power suppressed operators 
which are enhanced by a factor A 1 / 3 arising from the length of a large nucleus with mass number A [l7l. [l8l]. 

In semi-inclusive processes such as single hadron inclusive events in DIS on a large nucleus, a hard jet is formed in 
the collision of the virtual photon with a hard quark. This jet then begins to shower and lose virtuality on its way to 
hadronization. Some part of this space-time evolution occurs within the nuclear medium. Multiple scatteringof the 
jet in the medium modifies the final distribution of high momentum hadrons emanating from such a hard jet [19l,l20l]. 
Experimental measurements of single and multi-hadron production from such modified jets and their comparison with 
jets produced in DIS on a proton or in p — p collisions allow one to quantify the gluon distribution in dense extended 
QCD media 

This modification depends on a class of higher twist operators evaluated in the nuclear medium. While there exists 
considerable information regarding the ground state (nucleon) structure of large nuclei which may be invoked in the 
modeling of these higher twist operators, there exists practically no such information regarding the bulk structure of 
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the deconfined matter produced in high-energy heavy-ion collisions. The modification of hard jets, in the deconfincd 
matter produced in heavy-ion collisions, has assumed center stage in the experimental program as the primary probe of 
the structure of the produced matter [22| . This is due in part to the dramatically large effects seen in comparison with 
cold confined matter as well as the possibility for a first principles computation of this modification from perturbative 
QCD (pQCD). There are many approaches to this calculation, all involving a different set of approximations about 
the medium [23|, 0, US US H3] ■ For a review of the different approaches see Ref. [28[ ■ 

While the benefits of an effective field theory description of power corrections to hard processes in QCD cannot be 
overstated, to date there has not been a single attempt to incorporate the leading effects of medium enhanced higher 
twist within an effective theory formalism such as SCET. There exists an effective theory description of the dense 
deconfined matter in the limit of very high temperature. Here, a consistent effective theory of the medium without a 
jet is set up first, then the hard jet is assumed to have interactions with the soft field in the medium which are similar 
to those encountered by a hard thermal parton (29j . This article will take a different route, by trying to extend an 
existing effective theory of jets in vacuum to incorporate the effects of scattering in a medium. In what follows, we 
undertake the simplest extension to an SCET like formalism in the presence of an extended QCD medium. 

In Sec. II, the emergent scales in the problem will be discussed and the presence of a new mode, called the 
Glauber mode, not present in the current vacuum implementation of SCET will be motivated. In Sec. Ill, an effective 
Lagrangian which includes the interaction of such Glauber modes with collinear quarks will be derived from the 
QCD Lagrangian. In this first attempt, we will ignore the further interactions between these Glauber modes and 
the soft and collinear gluon modes of the usual SCET Lagrangian. A new set of Feynman rules arising from such 
an effective Lagrangian will be outlined and their equivalence with the Feynman rules of full QCD demonstrated at 
an amplitude by amplitude level. In Sec. IV and V, the Feynman rules will be used to compute cross sections in 
physical processes; two examples will be dealt with: the transverse broadening of jets in DIS on large nuclei and 
the transverse- momcntum-dependent-parton-distribution- function (TMDPDF). It will be shown explicitly that, in 
light-cone gauge, the Glauber gluons give rise to the transverse gauge link that enters the definition of the fully gauge 
invariant TMDPDF. The results obtained will be compared with published calculations in full QCD. Concluding 
discussions will be presented in Sec. VI. 



II. THE ENERGY SCALES FROM A DENSE MEDIUM 



Consider the deep inelastic scattering (DIS) of a hard photon with momentum q and virtuality q 2 = —Q 2 on a 
nucleon with momentum p in vacuum or contained within a large nucleus. In the Brcit frame, the photon has the 
momentum components 

q=[q + ,q-,qx] = [Q 2 /2<T, <T, 0] ~ Q(-l, 1, 0). (1) 

The off-shellncss of the photon Q is taken as a representative of the hard scale in the process. This strikes a hard, 
almost on-shell quark with a large momentum in the +z direction or a large (+)-component of momentum, 

Pi = {x B p + ,p-,p t ,±) ~ 0(1, A 2 , A). (2) 

The number of partons in the infinite momentum frame which carry an xb ~ 1 fraction of the nucleon momentum is 
rather small and A is a small dimensionless variable (A — > 0). We then trigger on events where an almost on-shell jet 
is produced in the final state (see Fig. [I}. This partonic jet moves with large momentum in the (— )-direction, 

p f = (p+,q-,p ftl _)~Q(\ 2 ,l,\), (3) 

where, at leading order, pj\±_ — Pi.j_- 

In the case of DIS on a single nucleon, this partonic jet immediately escapes the medium and eventually after a 
time ~ 1/(X 2 Q) 1 would have decayed into multiple partons of lower invariant mass and eventually turns into a jet of 
hadrons. The invariant mass of the final jet is m ~ XQ and its total forward energy from all produced hadrons which 
have arisen from this jet is 

E = q~/V2~Q. (4) 



The formation time of the radiation may be estimated from the virtuality of the jet ~ AQ and the boost of the jet ~ A . 




FIG. 1: Single-inclusive DIS on a nucleon or a large nucleus. 



In this article, we initiate the construction of the effective theory for the propagation of such jets through the dense 
matter within nucleons and large nuclei. 

The first step in such an endeavor is the assignment of relations between all relevant dimcnsionful quantities that 
appear in the problem. We assume that the scaling variable A is so chosen that perturbation theory may be applied 
down to momentum transfer scales at or above X 2 Q. The medium, also introduces its own set of scales, e.g., the mass 
of a nucleon Mjv ~ lGeV. This is assumed to scale with a new scaling variable \i such that it is comparable to Aqcd 
and in general much smaller than the soft perturbative scale, i.e., 

M N ~ A QCD ~ M Q < X 2 Q. (5) 

One may immediately surmise that in the Breit frame when ig ~ 1, and p + ~ Q, the boost or 7-factor is of order 
jjL . While the inverse size of the nucleon may be thought of as an even softer scale, in this effort, we will assume 
that the hard scales Q, XQ are much harder than the medium scale fiQ and thus the inverse length will be assumed 
to be of the order of the mass of the nucleon, 

In ~ 4v ( 6 ) 
fiQ 

All our considerations will be carried out in the Breit frame. The scaling introduced at the beginning of this section 
regarding the momentum components of the incoming and outgoing partons were set up in the Breit frame. Thus 
the nucleon (or medium) will have to be boosted to this frame. The ensuing boost to the Breit frame will lead to the 
contracted length of the nucleon, 

^ ~ -. (7) 

The very introduction of an alternative soft scale such as [iQ may lead the reader to imagine a much more complicated 
effective theory which will manifestly involve both A and fi and will require a relation between these two scaling 
variables. However, as will be demonstrated in the next two sections, with specific examples, it is possible to construct 
an effective theory in dense matter involving only the hard scale Q and the vacuum scaling variable A. 

The construction of an in-medium effective theory which depends on only the hard scale and the scaling variables 
from the vacuum theory has one further requirement. In specific examples, such as in Sec. IV, certain in-medium 
matrix elements will be enhanced by media with sizes much larger than a nucleon, e.g., in the case of large nuclei or 
a deconfined quark gluon plasma (QGP). In the case of large nuclei (with mass number A I), the enhancement 
factor is usually the nuclear length in units of the nucleon length i.e., A 1 / 3 . These situations will require that the 
enhancement be expressed in powers of A, i.e., A 1 / 3 ~ A - ™. The number n is so far unspecified and will turn out to 
be observable dependent. 

The particular choice of scaling of nucleon size and momentum lead to certain obvious physical consequences: In 
the Breit frame, valence (large x) partons carrying order one fractions of the forward momentum of the nucleon, have 
momenta that may be expressed as (by simply boosting momenta of the order of M/v), 

fc~Q(f,A 2 ,A), (8) 



These partons will be found to be completely confined within nucleons. The off-shcllncss of these partons ~ X 2 Q 
is very small and, as a result, for processes which involve momentum transfers of the order of XQ or larger, these 



4 



radiated partons may be considered as asymptotic in-states. In the following, we will often ignore discussion of 
the (— )-components of the partonic momenta in the in-state; these play almost no role in the computation of jet 
modification and the transverse momentum dependent structure functions. 

Due to interactions, a variety of gluons may be radiated from these valence partons with momenta constrained by 
overall energy-momentum conservation and determined by the kinematics of the process being triggered on. In the 
case of transverse broadening, the type of radiated gluon which plays a leading role will be those with momenta which 
scale as 

k~ Q(A 2 ,A 2 ,A). (9) 

Gluons with momenta which scale as in the above equation are referred to as Glauber gluons, or gluons in the Glauber 
region. The role of Glauber gluons in transverse broadening (as well as in transverse momentum dependent structure 
functions) may be easily understood in the case of DIS with a hard jet in the final state. The momenta of the produced 
quark jet scales as in Eq. Imagine the multiple scattering of the struck quark off the remnants of the nucleon or 
nucleus. The diagrams under consideration are of the form of Fig. O In order for the produced jet to escape from 
the nucleon or nucleus without undergoing any induced radiation, the interactions with the nucleon have to be such 
that they do not induce a major change in the off-shcllncss of the quark. In order to see how this comes about we 
explicitly write out the expression for the qi propagator where qi = q + pi + ki, represents the quark propagator, 
which, after the hard scattering (with momentum q + Pi) scatters off one extra gluon with momentum k\, (in order 
to simplify the expression we assume that p\ = and p\ = xbP + = —q + ) in Fig. [2j 



S(q2 



7 + q~ + l ± k]_ +7" kf 
2q-(k+)-\kl\2 

7 + g~ 
2q-(k+)-\kl 
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(10) 



Since, q~ ~ Q and k]_ ~ XQ, the forward momentum has to scale as ~ X 2 Q for the jet to remain off-shell by no 
more than \ 2 Q 2 . If the forward momentum scales with a higher power, i.e., k + ~ XQ, this will cause the jet to go 
off-shell by XQ 2 and lead to the radiation of momenta with large transverse momenta l± ~ A 5 Q. This process will 
lead to the radiative energy loss of the propagating quark and will be discussed in a future effort. The reader will 
note that the absorption of gluons which are collincar to the outgoing quark, i.e., with a momentum that scales as 
Q(X 2 , 1, A) also do not raise the off-shcllncss of the quark beyond X 2 Q 2 . However, the number of such gluons emitted 
from a medium moving with a large collinear momentum in the (+)-direction is vanishingly small. Hence the effect 
of such partons will be ignored. 



^2 
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FIG. 2: The multiple scattering of a produced jet in Deep-Inelastic scattering. 



In the current manuscript, the focus will remain on the production of a single jet with non-zero transverse momen- 
tum. Since the number of gluons with a forward momentum p + ~ X 2 Q far exceeds those with p + ~ XQ for the same 
transverse momentum, the jet will tend to encounter multiple interactions with gluons with a soft forward momentum. 
These will result in the transverse broadening of the hard jet. The neglect of gluons with a larger (+)-component of 
momentum suppresses radiation from the hard parton. 

With the power counting of the different momentum components identified, there remains the issue of determining 
the power counting (in terms of A) of the 4-vector potential A£ in this regime of momenta. In the case of effective 
theories of QCD in a vacuum, the power counting of the A% field is determined by an estimation of the powers of 
A from the gluon propagator. In the case of Glauber gluons such a methodology will yield incorrect results. In the 
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Glauber region of momenta, the gluon propagator obtained from the full QCD Lagrangian is never on-shcll.With the 
transverse momenta being larger than the light-cone components, Glauber gluons are always space-like off-shell. As 
demonstrated above, an shell collincar parton may interact with a Glauber gluon and have its transverse momentum 
changed by order XQ while still remaining on-shell. The simplest extension of an effective theory containing only 
collinear modes (the jet) to a medium with collincar modes travelling in the opposite direction (the target) will contain 
the interactions of the collinear jet parton with Glauber gluons radiated off the partons in the target which move in 
the opposite direction. The effective action has the simple form, 

S= Jd 4 x [L SCET +j^x)A a G (x)] . (11) 

where, A G is the Glauber field radiated from the target and is the current of the collinear partons from the jet. 
The kinetic and interaction terms for the collinear fields which constitute jfj are contained within Lscet 0] along 
with terms for the soft fields. In this effort, ^scet will not contain the collinear or soft modes from the target. These 
will be integrated out and included in effective Glauber field Aq. 

Note that there is no kinetic term for the Glauber gluon; thus, it does not obey a classical equation of motion. It 
admits no mode expansion and is not quantized as the SCET modes. However such exchanges are included in the 
full QCD Lagrangian and are prevalent in the interaction of collinear modes from jet with those from the target. 
The Glauber field A G represents the effective classical field of the target partons. The power counting of the various 
components of the Glauber field may be obtained from a calculation of its production in full QCD. 

As our goal here is simply to estimate the power counting of the various components of the Glauber field, we will 
ignore subtleties associated with non-linear terms in an interacting non-Abelian theory. We will estimate the A power 
of the various components in an Abelian theory. In a classical Abelian theory, the gauge field A^ is obtained from a 
solution of the inhomogeneous Maxwell's equation. This is given as, 

A»(x) = A%(x) + J d 4 yW(x-y)J u (y), (12) 

where, Aq(x) is a solution of the homogeneous Maxwell's equation. By restricting the current to be collinear to the 
target direction, and insisting that the incoming partons in the target remain close to on-shell, we restrict the field 
A fl (x) to only its Glauber component A G (x). In this region of momenta, there exists no solution of the homogeneous 
Maxwell's equation i.e., Aq g (x) = 0. As a result, the Glauber field is obtained from the second term on the right hand 
side of Eq. (fl"2"|) . We now evaluate, the power counting of the various components of the Glauber field in covariant 
and light cone gauge. 



A. Covariant Gauge 



At leading order in covariant gauge, the gauge propagator CD MI/ is given as, 

° A.. • (13) 

In Eq. (|12p . J u {y) = , 4 , (y)l v '4'(y) is the current of partons in the target which generates the gauge field. The fcrmionic 
operator may be decomposed as, 

iP(y) = f dp+d2p± - -. ]T u s (p)a;e-^y + v s (p)b^e^y (14) 

J (2^)3 - 

The scaling of the fermionic operator depends on the range of momentum which are selected from the in-state by 
the annihilation operator. Note that this influences both the scaling of the annihilation operator a p as well as the 
bispinor u(p). The power counting of the annihilation operator may be surmised from the standard anti-commutation 
relation, 

{a r p ,a s p S} = {2K) s 5 3 (?-if)8 rs . (15) 
and the power counting of the bispinor from the normalization condition, 

^2 uS p Up=P' = J~P + +J + P~-J±-P±- (16) 
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Substituting the equation for the current in Eq. (|13[) . and integrating out y, we obtain the expression for the (+)- 
component of the gauge field: 

i(p-q)-x 

a\a p u(q)"f + u{q). (17) 




If the incoming and out going momenta p and q scale as collinear momenta in the (+)-direction, i.e., p ~ Q(l, A 2 , A), 
then we get, 5 3 (p — p 1 ) ~ [A 2 <5 3 ] -1 , as one of the momenta will involve the large (+)-component and the remaining 
are the small transverse components. Thus the annihilation (and creation) operator scales as X~ 1 Q~~ 3 / 2 . Also in 
the spin sum f> ~ Q and thus u(p) ~ u(q) ~ Q 1 ! 2 . The 7 + projects out the large (~ Q) components in u and u 
in the expression u(q)j + u(p). We also institute the Glauber condition that p + — q + ~ A 2 Q, p~ — q~ ~ A 2 Q and 
P_l - gu_ ~ AQ. 

Using these scaling relations we correctly find that the bispinor scales as XQ 3 ^ 2 . However, to obtain the correct 
scaling of the gauge field A + one needs to institute the approximation that q + = p + + k + where k + ~ X 2 Q. This 
condition is introduced by insisting that the (+) momentum of the incoming and outgoing state, which control the 
scaling of and a p , are separated by k + ~ X 2 Q. This is used to shift the dq + — > <ifc + and as a result we obtain 
the scaling of the A + field as A 2 Q. Following a similar derivation, with the replacement 7+ — > 7 (7") we obtain the 
scaling of the transverse and (— )-component of the gauge field as A 1 - ~ X 3 Q and A" ~ A 4 Q. 



B. Lightcone Gauge 

The power counting of the gauge field is gauge dependent. In this last subsection we surmise the power counting, 
in light-cone gauge, for the Glauber field. The primary difference with Eq. (|17[) is the gauge field propagator. In the 
positive light cone gauge: n ■ A — n~~A + = A + = 0, the only non-zero components are A± and A~ . Note that a 
Glauber field with transverse momentum kj_ ~ AQ can only be radiated from a collinear parton without changing the 
direction of the collinear parton. For A±, the dominant contribution to the power counting equation arises not from 
the g^ v term in the numerator of the propagator, but rather from the {k^n v + k v n^)/k + term, i.e., 

. ,3 ,3 i( (P-L-gJ-)"" ^ e -i(p-q)-x 

w ^« 1 ' :: L — «i^hM«y us) 



(p - qf 

Comparing this with Eq. (fT7|) . we obtain that A± ~ AQ. Similarly we obtain A" ~ A 2 Q. As a result, in light cone 
gauge, the (_L)-component of the gauge field is much more dominant than in covariant gauge and we expect this to 
change the power counting of various terms in the effective Lagrangian. Similar power counting arguments may also 
be surmised from the explicit expressions presented in Ref. [39( and references therein. See also the discussion at the 
end of section V. 

In the next section, the power counting arguments presented above will be used to derive the effective Lagrangian 
which describes the interaction of collinear modes with Glauber exchanges with the medium. While these power 
counting arguments have been derived for an abelian theory, we expect them to remain true in a non-Abelian theory 
as well. 



III. EFFECTIVE LAGRANGIAN FOR GLAUBER GLUONS 



In the preceding section, the momentum scales associated with a new mode which arises in the presence of a 
medium was outlined. These, so called Glauber gluons present a mode that was absent in the derivation of the SCET 
Lagrangian. In what follows, we introduce these modes and construct a new additional effective Lagrangian called the 
Glauber Lagrangian. In this first attempt, the kinetic terms which represent the soft and collinear gluons of the SCET 
Lagrangian, along with their interactions with the collinear quarks will be ignored. The two active fields will be the 
collinear quarks and the Glauber gluons. In principle there will be a similar contribution from Glauber interactions 
with a collinear gluon. While such interactions are not included in the derivations presented in the current paper, 
these represent a simple extension of the formalism presented in this section. 

Before we start the derivation of the effective Lagrangian we comment on the off-shellness of the Glauber gluons. 
Since for Glauber gluons the product k + k~ is much less than |fcjj 2 then the Glauber modes arc obviously off-shell 
degrees of freedom. In principle when one constructs and effective Lagrangian, the degrees of freedom involved have 
to be on-shell so that one can make use of the classical equations of motions. In our derivation below we do not make 
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any use of the gluon equation of motion and it is only the Dirac equation for a collincar quark that is utilized. Thus 
the derivation below should be viewed merely as the limit of the contribution from gluons with arbitrary momentum 
taken to the Glauber region. This is much similar to the well-known method of regions. In this sense we are deriving 
actually an effective "vertex" between a collinear quark and a Glauber Gluon. This effective vertex could also be 
obtained on a case by case basis starting from the full QCD amplitudes and then taking the Glauber limit. This is 
illustrated below for a non-trivial case. The advantage of deriving an effective Lagrangian is mainly the consistent 
power counting (of the gluon fields and momenta) invoked in the derivation. 

In DIS, in the Breit frame, the final state after the hard scattering consists of a hard out-going quark in the — z 
or simply the (— )-dircction, which interacts with the remnants of the proton (or nucleus) which move in the +z or 
simply the (+)-dircction. These interactions are dominated by soft grazing scatterings with the fast moving remnants. 
Some of these interactions may be hard as well, resulting in the hard quark going considerably off-shell and radiating 
a hard gluon. Such processes constitute the energy loss of the hard quark and will be ignored in this first attempt. 
The incorporation of such interactions will necessarily involve the reintroduction of both the soft and collinear modes 
as well as the inclusion of new interaction terms between the Glauber modes and these soft and collinear modes. In 
what follows we focus solely on the soft interactions of the hard outgoing quark with the soft glue field generated by 
the remnants of the struck proton (or nucleus). 

Consider a fast-moving quark moving along the n direction where n = ^(1,0,0,-1). The +z direction will be 

denoted by the unit vector, n = ^75(1, 0, 0, 1). This quark has large momentum in the — z-dircction p~ = n ■ p ~ Q 
where Q is the hard scale of the process considered. The interaction of this collinear quark with Glauber gluons leads 
to change of the transverse momentum component of the collinear quark while the p~ component remains fixed up to 
0(A 2 ). Thus the truly label-changing component is only the transverse one p±. As such the starting point to describe 
the interaction of this collinear quark field with Glauber gluons would be, as in SCET @, §, to extract the label 
momentum components from the full QCD field. 

iKar) = e-^ + ^e^-^n,i>x(z)- (19) 

Again as in SCET, we decompose the "4>n,p ± hito a sum of two fields: £,n,p ± + £,n,p± where £n,p± carries the large 
momentum components in the — z direction while t; n ,p ± carries the small momentum components. The next step is to 
substitute this decomposition into Eq. ([19)) and then substitute the result into the interaction term in the full QCD 
Lagrangian. In order to maintain consistent power counting in the effective theory one has to specify the scalings 
(in terms of A) of the relevant quantum fields. For collincar quark the scaling is the same as in SCET, namely S,n,p ± 
scales as A while £ n ,p ± scales as A 2 . All derivatives acting on £n,p ± or £, n .p ± will further suppress the power counting 
by A 2 . For the Glauber gluon gauge field the scaling was given in the previous section. 

The starting point to obtain the effective Lagrangian is, as in SCET, the full QCD quark sector expressed in terms 
of the fields and £ n 

£g = £h,Px fti in ' D )£fi,P± + €n,P± f>-( n ' P+ in ' D )£n t f x 

+L,P' ± + 1 + L,p\ (*1 + * PlW ■ ( 2 °) 

We notice that in the last result there are terms that scale as 0(A 4 ), 0(A 5 ) and 0(A 6 ). We now eliminate the 
non-dynamical field £„ by making use of the tree-level equation of motion 

fi± + i pj_ 
2(n - p + in - D) 

We again have kept leading and sub- leading contributions in the Eq. (f2~T| . It is useful to notice the difference between 
the case of Glauber gluons and collinear gluons. For collinear gluons, the gauge field component in the covariant 
derivative (n ■ A) in the denominator of Eq. (f2~Tj) scales the same as n ■ p. This component of the gauge field, 
eventually, leads to the presence of the collinear Wilson lines in SCET, as was demonstrated in Ref. [{|. However, for 
Glauber gluons this covariant derivative is suppressed compared to n ■ p. Therefore, we expand the denominator and 
get 



€™,P± = K7Z Z I 7Z n\ ^n,p± ■ ( 21 ) 



£ 9 = E e ^~ A) ^w<>) 

p±,p'± 



n-iD + (j ± +i P^t-^— (l - ^-^ ) (jrfj. + i p±) 
In ■ p \ n ■ p 



fcnjAx) ■ (22) 



where higher orders in A have been dropped out. In this article we will only consider Glauber gluons in covariant and 
light-cone (A + = 0) gauge. In covariant Gauge the leading order Lagrangian is given by 

[ _ p2 -1 

£g = S,n,p ± in- D + fi£n,p ± (23) 
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where both terms in the square brackets are of order A 2 . In light-cone gauge the leading order interaction Lagrangian 
is given by 



2n ■ p 



(24) 



The Fcynman rules derived from the effective Lagrangians above are given in Fig. [3l As a first simple test of these 
rules we compute the amplitude for the case of two Glauber gluons attached to a collinear quark line and compare 
with the amplitude obtained from full QCD. We show that the effective theory exactly reproduces the full QCD result 
at the level of the amplitudes of Feynman diagrams. 



n-p 
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FIG. 3: Feynman rules for a collinear quark interacting with Glauber gluons. In the second diagram with a single gluon vertex, 
the first term is the contribution from the covariant gauge Lagrangian while the second term is the leading contribution from 
the light-cone gauge Lagrangian. The third diagram is only required for calculations in light-cone gauge. 

Let us consider the Feynman diagram given in Fig. [4] where two Glauber gluons are attached to the collinear quark 
field. In full QCD, the amplitude reads 

T--f ^2 f d 4 k, d 4 k 2 Ajk^p 1 - h\ h- fej / 9 -x 

[Ws) J (2tt)* (2^ U[P) [(p - k,r + ie][(p -h- k 2 f + is) ' [Zi)) 

where u(p) represents the Dirac spinor for an outgoing quark and p has no transverse momentum. The scaling of 
the quark spinor will be ignored in the following, as it plays no role in the remaining discussion. When expanding 
the numerator in Eq. (|25[) one should invoke the same power counting for the gluon gauge fields and the momenta as 
the one used in deriving the effective theory. In light-cone gauge A + = and by making use of the Dirac equation 
[u{p) p 1 = p~u(p)j + = 0], the leading contribution is, 

2 r d 4 h d 4 k 2 Axjkj^kfp- A±(k2)+hi- AJMKhx± h±)\ l0R , 

3 = - {l9s) J J2^J2^ U{P) [(p-k 1 y +l e][(p-k 1 -k 2 r +l e] ' (26) 

We notice that each contribution in the square bracket scales as A and subleading terms have been dropped out. 

In the effective theory (and again working in light-cone gauge) there are two Feynman diagrams that contribute. 
One (denoted below as 3^) comes from two Glauber gluons attached at different points. The other contribution 
comes from the vertex of two glauber gluons attached at the same point. It is denoted by 3^- Using the Feynman 
rules given in Fig. [3] we obtain the first contribution from the effective theory as, 

(1) = 2 r d% - AjM h±[A±(k 2 )(h±+ hx 4±(fc 2 )] (27) 

[Ws) J (2tt)4 (2tt) 4 U [ 2 p-k+ + |fc 1± P - ie][2p-(kt + k+) + \k ix + k 2 i_\ 2 - ie] ' 
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FIG. 4: Final state interactions in DIS: two gluon exchange. 



The second contribution from the effective theory is given as, 



a(2) ( \* f d ' kl rf4fc2 p A±(k ^ 

d ~ [lgs) J (fcr)* (2*)« ^ [ 2p -(fc+ + fc+) + \k lx + fc 2 J2 _ l£] ■ W 

fn the above result we have considered only one contribution where the color and Lorentz indices are held fixed in 
the Feynman rule for the two-gluon vertex. The other contribution gives an identical result as 3^ upon integrations 
over k\ and k 2 - It is can be easily verified that # = thus confirming that the effective theory reproduces the 

full QCD result at the level of the amplitudes of the relevant Feynman diagrams. The case of one-gluon attachment 
is trivial and one can easily verify that the effective Lagrangian also gives the same result as the one in full QCD. 

In what follows, we investigate two physical applications of the derived effective Lagrangian: The transverse broad- 
ening experienced by hard jets in DIS on a large nucleus and the transverse momentum dependent parton distribution 
function (TMDPDF) in a nucleon. In the first case, a final physical cross section will be computed and arguments on 
the enhancement of power corrections by large lengths in nuclei will be forwarded; hence, this application contains 
arguments beyond those used to derive the effective Lagrangian. However, the scaling of the momenta of the gluons 
will always lie within the strict boundary prescribed by the Glauber Lagrangian. It may not come as a surprise that 
the Glauber gluons which lead to the transverse broadening of hard jets, also play a principal role in the construction 
of the gauge invariant TMDPDF. 



IV. APPLICATION I: TRANSVERSE BROADENING IN LARGE NUCLEI 



A straightforward application of the effective Lagrangian derived in Sec. Ill, is to the process of jet broadening in 
dense matter. As a specific example, we consider the process of jet broadening in Deep-Inelastic scattering in large 
nuclei. A virtual photon with momentum q = [Q 2 /2n ■ q,n ■ q,0, 0] is incident on a large nucleus (A) with momentum 
Ap ■ n where p ■ n is the mean momentum of a nucleon. In the remaining section, we will refer to q ■ n as simply q~ 
and p ■ n as simply p + . 

We compute the cross section for the semi-inclusive production of a hard jet in the final state with a net transverse 
momentum lj_ with respect to the direction of the virtual photon, i.e., 

e(L 1 )+A(p)^e(L 2 ) + J(l ± )+X. (29) 

In the frame chosen, the Bjorken variable is defined as xb = Q 2 /(2p + <7~). The differential cross section may be 
decomposed into a leptonic and a hadronic part as, 

E L2 dcj a% M T dW^ u 



V-JoT-- (30) 



d 3 L 2 d 2 l 1 _ 2ttsQ 4 ^ d 2 l ± 

where s = (p + L\) 2 is the total invariant mass of the lepton nucleon system. The leptonic tensor may be expressed 
as, 

Lfiv = ^Tr[^ l7M £ 27 „]. (31) 
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The initial state of the incoming nucleus is defined as \A;p). The general final hadronic or partonic state is defined 
as \X). As a result, the semi-inclusive hadronic tensor may be defined as 

W^=Yi27r i )6 4 (q + P A -p x )(A;p\J^Q)\X)(X\r(0)\A;p) = 21m 
x 

where the sum (Ylx) runs over all possible hadronic states and J' 1 is the hadronic electromagnetic current i.e., 
= Qq^aj^^n, where Q q is the charge of a quark of flavor q in units of the positron charge e. It is understood 
that the factors of the electromagnetic coupling constant have already been extracted and included in Eq. (|30|) . The 
leptonic tensor will not be discussed further. The focus in the remaining shall lie exclusively on the hadronic tensor. 

In a full QCD calculation of Eq. |32|) . one computes the hadronic tensor, order by order, in the strong coupling. 
This leads to the introduction of a variety of processes leading to a modification of the structure of the jet. Such 
processes include radiative branchings, flavor changes of propagating partons, as well as transverse diffusion of the 
partons in the shower which ensues from the quark produced in the hard scattering. In this article, we will focus 
solely on the processes which lead to the transverse momentum diffusion or transverse broadening of the produced 
hard quark. 

In Ref. [30l ] . the leading contributions to transverse broadening without induced radiation, at all orders in coupling, 
were identified as those of Fig. [5j These diagrams depict processes where the propagating parton engenders multiple 
scattering off the glue field inside the various nucleons through which it propagates. However, scatterings do not 
change the small off-shellness of the propagating parton; as a result, large transverse momentum radiations do not 
occur. Using simple kinematics, the relation between the momentum components of the glue field fc,; may be surmised 
by insisting that the off-shellness of the i + 1 th quark line be of the same order as the i th line, 

(p + hf =p 2 + kf + 2p + kr + 2p-k+ - 2p x ■ k\. (33) 

Insisting that (p + ki) 2 ~ p 2 ~ X 2 Q 2 and given the known scaling of the quark momenta (i.e., p + ~ X 2 Q,p~~ ~ Q,p± ~ 
XQ), we obtain that k]_ ~ XQ, kf ~ X 2 Q and fcr may scale with a range of different choices Q, XQ, X 2 Q etc. The first 
two cases for the scaling of k~ represent gluons which are emanated with large (— )-momentum from a nucleon moving 
with large (+)-momentum. The number of such gluons must be vanishingly small. The first non-trivial population of 
gluons emanating from a nucleon moving with a large (+)-momentum, are those which scale as k ~ [A 2 , A 2 , A], which 
essentially constitute the Glauber sector. 



(32) 



q 



q; q 2 ' q,' 



Po l\ 



AP 




y/ y/ 




y 3 yi y y 



FIG. 5: An order n diagram which contributes solely to transverse broadening. 
Using the Feynman rules derived for Glauber gluons in section 2, the leading component of n th order diagrams such 



11 



as those of Fig. \5\ may be expressed as 



n 



n+n 



2q~k 



\kl\ 2 



1 



2l-2n5{2l+l- - J 



x t a M+ (^) 7 -|A;p) e -. 



j=n' 

let 



2g" 



2q-k' i + - \9 J x \ 2 +ie 



t ai A+( yi ) 



(34) 



where, it is understood that y i s the origin and y° = y. In the equation above, the gauge fields have been expressed 
in coordinate space. At this point an n th momentum may be introduced, via 



1 = J d 4 k n d 4 ^-X>fc-?J 



(35) 



This leads to a considerable simplification of the phase factors. The complete absence of the (— ) components of 
the momentum, from all expressions except for the phase factors allows for the k~ and k' integrations to be done, 
resulting in the localization of the process on the negative light-cone. 

The integrals over the momenta kf , fc' • may be re-expressed in terms of momentum fractions, i.e., 



Q = 2xbP q ; k l = Xip + ; k j = x } p 

i 

Y,2k i L -ki + \ki\ 2 = 2x i D p+q- 

k=0 
3 . , 

^2fc'i • k' ± + \k'i\ 2 = 2x' 3 D p+q- 

1=0 



(36) 
(37) 

(38) 



Integrating over all the Xi and x'j momentum fractions by contour integration, we obtain the much simplified form of 
the hadronic tensor, 



= g n+n ' 



v ; i=0 1 = 1 J i=0 K ' 1=0 v ' 



(2tt)V(Z 



n^r-%-i)II%7- 2/7-i) 



i=0 



IF 

3=0 



3=n 



x (A; P \Uy-,y±h + U0)Tr 



i=i 



\A; P ). 



(39) 



The expression derived above has so far been a direct application of the Feynman rules derived in the preceding 
section. Hitherto, no assumption regarding the nature of the nuclear state has been made. As a result the nuclear or 
nucleon scale of [iQ has also not appeared in any of the expressions. However, the hadronic tensor in Eq. (|39[) and 
any resulting transverse broadening will, ultimately, depend on the expectation of the (n + n' + 2)-parton operator as 
indicated in the last line of Eq. (|39[) . To proceed further, approximations regarding the expectation of this partonic 
operator will have to be made. In these approximations, the in-medium scale fiQ will appear. However, as we will 
show, the final transverse broadening will turn out to be independent of this scale under certain assumptions. 

Following standard treatments, we approximate the nucleus as a weakly interacting homogeneous gas of nucleons. 
Such an approximation is only sensible at very high energy, where, due to time dilation, the nucleons appear to 
travel in straight lines almost independent of each other over the interval of the interaction of the hard probe. All 
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forms of correlators between nucleons (spin, momentum, etc.) are assumed to be rather suppressed. As a result, the 
expectation value of the n 



A;p 



nf + 2 operators in the nuclear state may be decomposed as 

n+n 



A;p)=q 



P0,P2, 



(ri+ri')/2 

Pn {poMy-,y±-h + ZM \ ( Pi \A+( yi )Ai'.(y'i)\Pi)m 



i=l 



where, the factor 



represents the correlations between the (n + n')/2 nucleons which interact with the 



propagating parton. In the decomposition above, we have restricted at most two parton operators per nucleon, 
insistin g t hat any larger number of operators is suppressed. Note that this is only true outside the saturation 
regime j3ll . 32]. The decomposition performed above, also restricts n = n! . The choice of gluon operators per nucleon 
is also special to the case of transverse broadening: The maximum broadening is obtained when one gluon operator 
from the amplitude is paired with one from the complex conjugate. This reason for this is immediately understood 
with further simplifications on each gluon pair (written with spin and color indices suppressed), 



^ldYi(^ + (ylM + (2/l)b) cx P -ix* D p + y7 + ■ y± + '>>V, -ik'l ■ y'l 



(2w) 2 6 2 (kl - fc'i) / d 2 y ± e-^P + ^-y'7) e ^-y^(p\A+(y ± /2)A + (-y ± /2)\p), 



(41) 



where, y± is the transverse gap between the two gluon insertions and kj_ = (k^ + k' ± )/2. The physics of the above 
equation is essentially the transverse translation symmetry of the two gluon correlator in a very large nucleus. This 
is then used to equate the transverse momenta emanating from the two gluon insertions. Thus, if the two operators 
were both chosen from the amplitude or the complex conjugate, then the momenta brought in by one gluon operator 
would be immediately taken out by the other and, as a result, the combination of the two operators will lead to no 
net transverse broadening. The integration above, also simplifies the longitudinal phase factors which now depend 
solely on the difference of the longitudinal positions of the two gluon insertions. 

Further simplifications are introduced by Taylor expanding the transverse momentum dependent delta function, as 



n 



d 2 



\ 2\d 2 k\ 



5 2 {h 



n 



(42) 



and combining the |fcjj 2 with the expectation of the two gluon operator (pi\A^.(y l ± /2)A^,(—y l ± /2)\p i ) to convert 
these into the expectation of field strengths in the nucleon (p i \F+- L (y]_/2)F£- L (— y\/2)\pi). The meaning of this 
decomposition of the transverse momentum delta function is the retention of solely the leading twist part of each of 
the two-point correlators in each nucleon. Higher powers of a given transverse momentum will necessarily lead to 
higher transverse moments of the two gluon operator. One further assumption regarding the two point function of 
Eq. (|41[) , due to color confinement, leads to a constraint on the two longitudinal y~ integrations (ignoring color and 
spin indices), 



dy-dy' (p\F(y-)F(y' )\p) ~ / dy~(FF)y 



Icon ft 



(43) 



where, (FF) is the gluon expectation in a nucleon and y con f is the confining distance. 

Each such integral yields a factor of L~ ~ A 1 / 3 ~ 1/A from the unconstrained y~ integration. The equating 
of the pairs of transverse momenta that appear in each two-gluon correlation, as well as the relation between the 
longitudinal momenta from the 9- functions in Eq. (|39[) , require that the largest transverse momentum broadening and 
largest length enhancement arises from the terms where the gluon correlations are built up in a mirror symmetric 
fashion, i.e., where the gluon insertion at y 1 is contracted with that at y' 1 . 

Averaging over the spins and colors of the two point functions in each nucleon, the remaining n longitudinal position 
integrals for the gluon insertions may be simplified as 



Y[dvi0(vr -%-x) = rf / LK~ 

A — 1 ' J A — 1 



(44) 



Invoking the above simplifications, the leading length enhanced contribution at order 2n to the differential hadronic 
tensor is obtained as, 



d 2 W^ 
d 2 U 



= C, 



{v^rs 2 (Zx) 



2Nr. 



(45) 
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where Wq is the leading order transverse momentum integrated hadronic tensor, given as 

Wr = 2n[g^-g v + + g v ~ - <T] £ Q 2 q f q (x B ), (46) 



where, the expectation of the two quark operator in Eq. (|40|) . leads to the quark structure function in the equation 
above. 

There remains the overall coefficient C po Pn which contains the weak correlations between the various struck 
nucleons. A study of such correlations in Refs. [H, [34| revealed that a simple factorized form such as C pa p = 
Cp g (p/2p + ) n , where p is the density of nucleons in a nucleus, is not completely inappropriate. Using this simple form 
one may sum over all n, i.e., over multiple scatterings of the quark in the nucleus, to obtain the rcsummcd equation, 



d 2 w^ _ {DL - )Vl <fwr 



d 2 l ± d 2 l ± ' 

where, d 2 Wq" /d 2 l± = Wq" 5 2 {l^) , and the constant D is given as, 



^ ;v ii _u (47) 



n 2 a s f d 3 yd 2 k± 

~2N~ C P } (2tt) 3 2 P + «*. ™ — 1 " I 2 g- 



D = -^rPi CT^T (p\F a+a (y)F\+(0)\p) ex P { -i ( l^j-y- - fc^ ■ y ± ) }. (48) 



It is this constant D which controls the broadening experienced by the hard jet in the extended nucleus. 

As shown in Rcf. (30j . Eq. (|4"T|) is a solution of the two dimensional transverse momentum diffusion equation, where 
the initial condition may be taken as a 5-function in transverse momentum. Taking moments of the solution of the 
diffusion equation, we obtain the total transverse momentum squared acquired by the hard quark after traversing a 
length L~ in the nucleus as given by the simple relation, 

(k 2 ± ) L - = ADL-. (49) 

The reader will note that we have used a two dimensional delta function as the input to the diffusion equation. This is 
an approximation to a very peaked distribution and one may use any other input distribution as well. The net extra 
broadening experienced by the input distribution is given by Eq. ()49|) . Given that the initial parton is an SCET mode, 
the transverse momentum is of the order of |fcjj 2 ~ (XQ) 2 . As a last step, we will demonstrate that the broadening 
obtained from multiple scattering in the large nucleus is of this order in power counting and thus one may continue 
to think of an SCET mode propagating in the extended medium. 

The power counting of net transverse momentum squared may be easily estimated from counting powers of A and 
H in the expression for D in Eq. f|48[) . As the expression is frame independent it will be evaluated in the Brcit frame. 
In this frame, the z-component of the length of the nucleon (~ l/(fiQ)) is contracted to a length of order 1/Q; hence 
the nucleon density p scales as 

p = I ~ p 2 Q 3 . (50) 

The dimension of the nucleon ket is obtained from the standard normalization of the on-shcll nucleon state, given as, 

(p\ q ) = (27r)->+%+ - q + )5 2 (f ± - gl) \p) ~ (fiQy 2 . (51) 

The F +J -F +1 - correlator scales as (X 3 Q 2 ) 2 from the standard Glauber scaling rules for the transverse momentum and 
the vector-potential. The enhanced length in the nucleus may be expressed in terms of the nuclear parameter A 1 / 3 
and the length of a nucleon In, as 

L-=Ah N ~^. (52) 

Substituting the above relations in Eq. (|48[) . and noting that the y~ and the y± coordinates are conjugate to the 
\k±\ 2 /2q~ and p± momenta, yields the A power counting of the net transverse momentum squared picked up by the 
hard parton as, 

(p±)l- ~ A*A 4 Q 2 , (53) 

independent of the medium scaling parameter p. The broadening is rather small in an object the size of a nucleon, 
but may get enhanced in large or dense media. As a result, for small nuclei where A 1 / 3 <C A~ 2 one may ignore this 
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extra effect of final state multiple scattering. For nuclei, where A 1 ! 3 ~ A -2 , or the medium has a very large gluon 
density we may obtain a broadening which is comparable to the jets inherent transverse momentum. In this case, we 
are in the "SCET-Glauber" region, where the derived effective Lagrangian in this paper may be used in combination 
with the SCET Lagrangian to understand the interaction of hard jets in dense media. For nuclei where A 1 / 3 3> A~ 2 , 
the Glauber modes will broaden the propagating jets beyond the scaling assumed in the derivation of the SCET 
Lagrangian and a different set of effective theories will need to be constructed. 

There is an unknown quantity that has been invoked a number of times in the discussion above: the inherent gluon 
density. This is the density of "small x" gluons that emanate from the current density in the medium and interact 
with the hard jet. The number of such gluons is a dimensionless quantity and thus difficult to estimate in a power 
counting calculation. The number of such gluons may in general also depend on the media in question. It is well 
known that the number of such gluons may become rather large at high energies and may thus lead to considerable 
broadening of hard jets. While the application in this section has focused on the broadening of jets in nuclei, the 
factorization properties afforded by SCET bode well for the applicability of this theory to jet broadening in dense 
Quark-Gluon-Plasmas (QGPs) created in heavy-ion collisions. While QGPs have been estimated to be from 10 to 
100 times denser than nuclear matter, their lifetimes are rather short and the majority of jets propagate rather short 
distances in the densest part of such environments. Given these experimental considerations, we expect the derived 
effective theory in combination with SCET to have wide applicability. 



V. APPLICATION II: TMDPDF 



Inclusive hard scattering processes like DIS can be factorized, into perturbatively calculable short-distance quantities 
convoluted with non-perturbative long-distance quantities [l[. The latter quantities are the familiar Feynman PDFs. 
For semi-inclusive processes where a single hadron is observed in the final state with a given transverse momentum 
then it is the TMDPDF that enters into the factorization formula for the cross-section. More details can be found in 
[Hi!]. The TMDPDF have been introduced long time ago in 0], as, 

f(x,k ± ) = l -J ^ ) t ^~ i(rk+ ~^' i±) x TO(r,a)4 x (°°>£~) x 7 + io(oo, 0)V(0, 0±)\P), (54) 
where L is path-ordered gauge link along the light-cone in the (— ) direction, i.e., 

L f (oo,t-)=Pexp(ig B [ dCA + (C,Z±))- (55) 



The Wilson line has a well-known origin: It comes from radiation of gluons which are collincar to the incoming 
parton. In SCET these Wilson lines are also the familiar collinear Wilson lines W In this work we will not discuss 
further the emergence of those Wilson lines as they are not related to Glauber gluons. 

The TMDPDF is a physical quantity and thus has to be gauge-invariant under arbitrary gauge transformation. 
However the above definition is gauge-invariant only in the set of non-singular gauges like covariant gauges where the 
gluon field vanishes at £~ = do. In singular gauges like light-cone gauge (with .4+ = 0) the gluon field (specifically 
the transverse components A±) does not vanish at £ _ = oo and a gauge transformation performed with Aj_(£~ = oo) 
will generate a non-vanishing phase that is not compensated by any gauge link. Thus the above definition of the 
TMDPDF has to be modified by introducing an additional gauge link formed from the transverse components Aj_: 



L f -= oo (a,0x)=Pex P fi«7.y d£i-A x (r = oo,fx)J . (56) 

where the line integral in the transverse plane can be performed in arbitrary direction. The above observations were 
first made in Ref . [39( . 

The important question that arises is what kind of interactions, say in DIS, build up this gauge link. This question 
was answered in the work of Belitsky, Ji and Yuan (B JY) [40] . There it was shown that the final-state interactions 
between the struck quark and the remnants of the incoming proton are responsible for the appearance of this gauge 
link. Those final-state interactions are mediated by Coulomb gluons that carry mainly a momentum in the transverse 
direction. In the next section we will check whether these gluons are Glauber gluons or not. It is important to verify 
this issue as the final state interactions are responsible for many physical effects like single-spin asymmetries and 
shadowing [4l], [4^, |43| . Any attempt to formulate an effective field theoretic approach to study such effects has to 
start from identifying the relevant momentum modes that mediate the interactions. 



15 



We start by briefly reviewing the work of BJY by considering the Fcynman diagram given in Fig. [6] A quark 
propagator with momentum p — k has a denominator: —2p~k + + k 2 where the quark has essentially large momentum 
in the — z-direction. The integral over the gluon momentum k gets contributions from a vanishing denominator. This 
could happen if (1) k is collinear to the outgoing quark. In this case both terms in the denominator scale as A 2 or (2) 
k is soft which means k 2 -c k + and k + scale as A 2 or (3) k has Glauber scaling. Assuming the gluon is emitted from a 
fast-moving nucleon in the +z direction then (1) is highly improbable. Case (2) does not lead to any transverse effects 
like transverse gauge link or transverse broadening as all componets of k scale similarly. Moreover soft gluons give rise 
to the familiar Eikonal soft Wilson lines. In SCET soft contributions can be handled by field re-definitions [l2| and 
they get factorized from the collinear sector. The remaining contribution comes from (3) where p~k + ~ \k±\ 2 — k 2 . 

In BJY, the exact scaling of the gluons is not specified. This issue will be addressed in some detail below. For now, 
we continue the review of their work. By making use of the Chisholm's representation, one obtains the following form 
for the propagator, 



1 



2p~k+ + k\ - is 



dTe~ lT{2p fc+ + fc i-^) ; (57) 



where the left hand side is obtained from the limit r = 0. With the above representation one is then able to carry out 
the integrations over k + for the amplitude of Fig. [5J One then gets the gluon field A^(x) (in the mixed coordinate- 
momentum representation) evaluated at A M (x~ = 2rp~ ,x + = 0, k±). Now take the scaling limit p~ — > oo before 
performing the r integrations. This sets the argument of the gluon field in the — light-cone direction at infinity and 
all the remaining r-dependence is now in the exponent. Then perform the r integration. This will result with a 
|fcjj 2 in the denominator. By repeating the above set of manipulations one gets the following result for a multi-gluon 
exchange, 

( P -,n\uo)\p) = (-lranigxuip-m^im [ ^% F^ 1 ¥ . ^m\ p )- m 



( 2 -) 2 'lE} =1 £±l a -fe 

In light-cone gauge it is Eq. (f58|) that gives the transverse gauge link. It is important to notice that the last result 
could be simply obtained by calculating, in full QCD, the amplitude for a Feynman diagram with arbitrary number 
of gluon attachments and then setting all the k + components of the gluon fields to wherever they show up. The 
last equation is only valid if one takes the scaling limit first which amounts to setting all the k + to 0. This procedure 
clearly violates the Glauber scaling as one has to maintain the relative scalings of k + and \k\ 2 /(2p~) in tact. 

The transverse gauge link that results from final-state interactions was also derived with somewhat different set 
of manipulations in [471 ] however the basic observation is still the same: it is gluons with vanishing k + that give 
rise to that gauge link. Another important issue related to the correct definition of the TMDPDF was raised in 
Refs. [UEll. There it was claimed that a soft factor, built up of soft Wilson lines, needs to be subtracted from the 
standard definition of the TMDPDF so as to get the desired features of the anomalous dimension of the TMDPDF. 
Soft factor subtractions were also discussed in the traditional literature of perturbative QCD (pQCD) [49l. IBl)!. l5ll. [52l| 
where this subtraction is aimed to avoid double-counting among mainly the collinear and soft contributions. In the 
effective field approach the double counting issue was treated within the "zero-bin" subtraction [53[ and later on a 
connection was made with the pQCD one [5J, [H, [56| . It is interesting to see wether the arguments of soft subtraction 
based on the anomalous dimension arguments are equivalent to the double counting issue. We will not discuss this 
issue further here and we leave it to a future work. 



A. One Glauber Exchange 

We start our analysis of the TMDPDF in light-cone gauge n ■ A = by considering the relatively simple case of only 
one Glauber gluon interacting with collinear quark. The kinematics are that of DIS: the incoming quark is collinear in 
the +z direction with p + large and of order of Q. This quark carries a longitudinal momentum fraction x and carries 
transverse momentum fcj_. The virtual photon moves in the — z direction with momentum q — (0, 0,0, — Q). Thus 
after the hard interaction the quark has momentum with p~ = Qj \[2 and transverse momentum k± . The Glauber 
gluon has also kj_ so the outgoing quark has only p~ . Following the notation of BJY we consider the contribution of 
Fig. [S]to the matrix element (p~ , N\j l/ (Q)\P). 

For Glauber gluons we now use the Feynman rules given in Fig. [3l We get 

- r 1 - f dk+ r d 2 k ± it ± - 

2p J 2tt J (2tt) 2 k+ + \k_iZ_ _ i£ 
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FIG. 6: Final state interactions in DIS: one gluon exchange. 



where we used vM = u. Wc now invoke the Fourier transform of A(k + , k±): 



A^xik+M) = / dx- / d 2 f ± i M± (.T-,^)e l ( fe+a; ~-^-^ (60) 



and substitute for A(k + , kj_) into Eq. (|59p and carry the k + integral by countour integration picking up the pole from 
k + = —-^4 — h is. The result is (from now on we drop the tilde on ^4), 



— I 



h = ig s In^z J dx~9{x~) J d 2 xxA^(x- ,xa_) J fixe '\ ^ (61) 

In order to carry out the integral over d 2 k± we complete the square in the exponent and shift the integration 
variable to k' ± = kj_ + £=xj_. The resulting integral proportional to k' ± vanishes by symmetry. The remaining d 2 k± 
integral is obtained as, 



00 d 2 k ± _j^ 2 *- 2 P - i / f°° „ y 2 P 



-oo(^) 2 x-(2nf\J_ 00 J r (2,)2 V 'V2 



2 



— I p 

2?r ar 



(62) 



Introducing the above mentioned simplifications in Eq. (|6ip . the result for 1\ reads, 



h = 9s — J dx~9(x~) J d 2 x± ^±A fl± (x^ , x±) f — J e 2*- . (63) 
Let us consider the x~ integral: 

/= J dx-9(x-)A„ x (x-,x x )(^j e~* P J- ±|2 (64) 

and perform the integration by parts. To do so we notice that in light-cone gauge the x~ dependence of A^± on x~ is 
just a 8(x~) [23, [4^,[4al (see Eq. (j6"8|) below.) Moreover the contribution from a highly oscillating phase (obtained from 
the lower limit of the integral) would give a vanishing contribution upon integration over the transverse coordinates. 
The result of the integral over x~ is then given as, 

2ip~~ 

I = — —77 x jVl(oo,xj_). (65) 

\xi_Y 

Substituting the above result in the expression for ii, we obtain, 

1 f -. 0± 

h = igsTT^ I rf2 ^-L i- ( 66 ) 
2tt J \xj_\ z 
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In what follows, we decompose, x± as x± = \x_i_\no where fig = (cos 9, sin 9). The integration measure is given as 
d x± = \xj_\d\xj_\dO. Substitution in the expression for Ii leads to the form, 



^ /"27T />00 



h = Ws^l^l 1 J d9 J d\xx.\{n e ) i A llX {o ,x A _), (67) 

with i = 1,2. Using the trivial relation: 7 AtJ ~7 t = ^([7 M_L ,7 l ] + {7* 1 >7*}) and the fact that in light-cone gauge the 
gluon field A^± at x~ — > oo is a pure gauge [40(, we obtain, 

A- 1 (a; - ->■ oo,xj_) = 9(x-)V<j)(r), (68) 
where r = |xj_| and </> is an arbitrary scalar function. Then we have 

A ± (x~ — ► oo,a;j_) = ^-fig, (69) 
ctr 

which shows that the gluon field Aj_ is directed in the radial direction. With this, it is straightforward to show that 
the contribution from the commutator of the 7 matrices vanishes by symmetry. Thus we obtain, 

h = ig s -z- \ d9 d\x±\(he)iAj- (00, x±) = ig s / d\x±\n ■ A(oo,x±), (70) 
27r J a Jo Jo 

where, the second equality in the equation above is derived from the use of Eq. (|69j) . 

The above analysis shows that it is indeed the Glauber gluons, arising from final state interactions, that build 
the transverse gauge link. In this sense, any (gauge invariant) effective field theory formulation, of the TMDPDF in 
particular (see e.g. Ref. (Hz!) or semi- inclusive hadronic processes in general, requires the introduction of the Glauber 
mode in addition to the soft and collincar modes. 

Two remarks arc in order. We first notice that the above treatment could also be carried out with a multiple of 
Glauber gluons attachments. The power counting of such Feynman diagrams would still be a leading one since those 
contributions arise from the leading order Lagrangian. The sum of all those contributions would give the transverse 
gauge link. Secondly we notice that the power counting of the Glauber gluon field, in light-cone gauge could be 
read-off from Eq. (|69[) . In Feynman gauge, the power counting of A + could also be read-of from explicit expressions 
(see, e.g., Eq. (14) in \^.) 



VI. CONCLUSIONS 



Effective theories now constitute a mainstay in the collection of theoretical methods used to apply perturbative QCD 
to phenomenological questions. The Soft-Collinear-Effective-Theory has been identified as a rigorous and systematic 
effective approach in the application to phenomena involving hard jets in vacuum. In this article, we have instituted 
the first extension of this "leading twist" effective theory to include power corrections from the medium. As a guide to 
understanding the effects of the medium on a jet, we have focused on a description of the rescattering encountered by 
a hard quark produced in deep- inelastic scattering on a nuclcon in vacuum or within a large nucleus. While most of 
the results derived in this manuscript are immediately applicable to quark jets propagating through confined media, 
these may be straightforwardly extended to gluon jets as well as to propagation in deconfined media. 

A jet in SCET is endowed with a very particular relation as regards the range of its different momentum components, 

lf = lp + ,p-,0±]~Q[\ 2 ,l,\], (71) 

where, we have specified the case for a jet moving in the (— )-direction with Q, a hard scale and A a small parameter. 
The virtuality of this jet allows it to resolve modes in the medium with transverse momentum ~ XQ. If the 
forward or (+)-momentum components of these in-medium modes scale as Q (or even as XQ), this will result in an 
intermediate parton with large off-shcllncss of the order of Q 2 (or AQ 2 ) and almost immediate hard radiation with 
large transverse momentum. Such interactions will, no doubt, change the large momentum label of the propagating 
SCET mode and will be dealt with in more detail in a future effort. If the forward momenta scale as X 2 Q, the 
off-shellness of the propagating mode remains within the scaling prescribed by SCET and as a result, the simplest 
cxtcntion to this vacuum theory is suggested in Sec. II: the interaction between hard collincar quarks (or gluons) with 
gluons in the medium which scale as 



k» ~Q[X 2 ,X 2 ,X}. 



(72) 
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Such gluons are referred to as Glauber gluons and in Sec. Ill we have constructed the effective Lagrangian which 
describes their interactions with the hard collinear modes. Although we have denoted the scaling of the fc _ momentum 
to be A 2 Q, it could indeed have any scaling k~ Q i.e., not be a hard collinear mode traveling in the (— )-direction. In 
the Breit frame with the medium moving with a large boost 7 ~ A~ 2 in the (+)-direction, such modes are energetically 
disfavored. 

Since Glauber modes carry a small fraction of the forward energy of the nucleon begin struck by the hard jet moving 
in the (— )-dircction, they are quite pervasive and thus the inclusion of such modes and their interactions with the 
hard collinear modes is rather important. Such interactions occur continuously on a hard jet propagating through a 
dense medium. Given their off-shcllness, SCET modes may traverse distances of the order of (A 2 Q) _1 before decay. 
When hard jets traverse large distances in dense matter, their total transverse momentum distribution is broadened. 
As a first application of the effective Glauber Lagrangian, this transverse broadening is derived for the case of DIS 
on a large nucleus in Sec. IV. Multiple interactions with Glauber modes may eventually lead to the generation of 
off-shellness or transverse momenta beyond the range of applicability of SCET and a completely different effective 
theory will have to be constructed. The transverse broadening as a function of the distance travelled, allows for an 
estimation of the range in size of media, within which, the effective theory will remain applicable. This is estimated 
in Sec. IV, with the aid of some phenomenological input. It is argued that the derived effective theory has a wide 
range of applicability which may easily encompass jet propagation in the cold confined matter in large nuclei to that 
in hot deconfincd matter created in high energy heavy-ion collisions. 

In Sec. V, as a second example of the role of Glauber gluons in hard processes, we have considered the treatment, 
in full QCD, of Belitsky, Ji and Yuan for the TMDPDF. In their calculation, they have shown that gluons with solely 
transverse momentum components build up a transverse gauge link which should be an integral part of the gauge 
invariant definition of the TMDPDF. Their analysis seemed to depend on taking the scaling limit first. In the current 
effort, we have demonstrated that the transverse gauge link may also be derived by keeping the Glauber scaling of 
gluon momenta between k + and \k±\ 2 /2p~ explicit throughout the calculation. The full link structure, in any gauge, 
has been shown to arise from a combination of collinear and Glauber gluons. 

As a final remark we address certain situations where the Glauber gluons do not contribute. A standard example 
would be DIS on a nucleon with its related physical quantities: The quark form factor and the PDF. It has been 
demonstrated that the only relevant modes that produce the infra-red behavior of QCD for DIS are the soft and 
collinear (see, e.g, [58| and references therein). Also for the factorization of the PDF itself (in the large x limit) 
similar arguments and conclusions have been given in Ref. [k|. The fact that Glauber gluons do not contribute to 
the factorization of DIS on a nucleon could, in principle, be shown (in the effective field theory approach) when one 
combines the soft, collinear and Glauber in one framework and then certain cancellations of the Glauber contributions 
should become manifest. This is a somewhat more involved topic, which we leave for a future effort. 

One may also compute the corrections to the single hadron inclusive cross section in DIS on a nucleon from the 
Glauber sector. In large nuclei, the produced jets tend to have a distribution in transverse momentum which is much 
wider than the case of DIS on a nucleon. It is well established that such corrections are not leading and are power 
suppressed, hence are unimportant in the case of DIS on a nucleon. This result is also consistent with the Glauber 
Lagrangian derived in the current manuscript. The magnitude of the correction from Glauber scattering may be 
estimated from Eq. ([55]) . by setting A = 1. It is clear that the (pj_) generated is suppressed by A 2 compared to 
that in a purely SCET process. Thus in the computation of the single hadron inclusive cross section from DIS on a 
nucleon, the contribution of the Glauber Lagrangian is power suppressed. In the case of the Drell-Yan process, the 
relevance of Glauber gluons is a more complicated issue. We will not discuss it further and refer the reader instead 
to Refs. @, [I| HI, HI- 

In future efforts, the interaction between the Glauber modes emanating from the medium and the soft and collinear 
gluons of the SCET Lagrangian will have to be derived. This will represent the first complete theoretical description 
of jets with off-shellness in the range of (XQ) 2 propagating through dense media. The setup of such a formalism will 
allow for the first systematic approach to such difficult problems such as factorization in hard jet production and 
modification in heavy-ion collisions. Embellishments of the Heavy-Quark-Effective- Theory (HQET) with Glauber 
modes will lead to more rigorous formulations of heavy-quark propagation in dense matter. Such extensions of SCET 
and HQET will lead to important advancements in our understanding of parton propagation and energy loss in dense 
matter and will no doubt play a leading role in the detailed theory and experimental comparison currently underway 
in DIS and heavy-ion collisions. 
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